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The above seminvariants may be converted into semin variants of (A) 
by means of equations (5), (6), (7). A comparison of these equations 
with (9) shows that the desired seminvariants of (A) are obtained simply 
by replacing in the above semi-canonical forms r iki by u^i, vm by v m and 
*"m by w m . 

A comparison of the seminvariants I^ rsl ^ of (A) with the corresponding 
seminvariants 5 for the case n = 2 shows the former to be independent. 
Moreover, the functional determinant of l[ rs() with respect to u ik i for 
•each value of I = 0, 1, 2. . .,m — 3 shows that l\ rsi) are independent 
among themselves and of the seminvariants I^ rst \ Equations (10) show 
that we have the proper number of solutions for the variables involved 
and that all other seminvariants of the complete system can be obtained 
by the differentiation of I (r s l) and l\ r s ° . We have therefore the following 
theorem : 

All seminvariants of (A) are functions of I {r s ^ (r = 0,1. .,« — l;r + s < 
n; t=l, 2, 3; t < s), l\ rst) (r = 0, I. .,n - 1; r + s < n; t = 1, 2; 
t < S ; I = 0,1,. .m - 3), and of the derivatives ofI (rsi) (t = 1,2) andl { { si) . 

1 Wilczynski, E. J., Projective Differential Geometry of Curves and Ruled Surfaces, 
Teubner, Leipzig, Chap. I. 

2 Wilczynski, E. J-, Ibid., Chap. II. 

3 Stouffer, E. B., London, Proc. Math. Soc. (Ser. 2), 15, 1916 (217-226). 

4 Green, G. M., Trans. Amer. Math, Soc, New York, 16, 1915 (1-12). 

6 Stouffer, E. B., London, Proc. Math. Soc. (Ser. 2), 17, 1919 (337-352). 



SOME NEW METHODS IN INTERIOR BALLISTICS* 
By Arthur Gordon Webster 
CiyARK University, Worcester, Mass. 
Read before the Academy, April 26, 1920. 

The principal problem of interior ballistics is, given a particular gun, 
a particular shot, and a particular kind of powder, to find, for a given 
load, the position and velocity of the shot, the mean pressure (and inci- 
dentally temperature) of the gases in the gun, and the fraction of the 
powder burned, all as functions of the time or of each other until the 
exit of the shot from the muzzle of the gun. In particular, we wish to 
know the muzzle-velocity of the shot, the maximum pressure to which 
the gun will be exposed, and the portion of the bore which will be exposed 
to it. It is then the duty of the mechanical engineer to design a gun to 
safely resist the pressure that may be expected. Or it is the inverse 
problem of the ballistician, by experiments on the action of the powder 
in the gun, to find its properties and those of the gun. It is true that the 
properties of the powder are more conveniently studied by means of 
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explosions in the bomb at constant volume, but we shall here show how 
all the properties may be determined by experiments in the gun alone. 

It is customary among ballisticians to make many more or less crude 
hypotheses, and it is even denied that it is possible to apply the strict 
principles of thermodynamics to the gun as a heat-engine. It shall be 
our endeavor to show that this is not the case, and to relieve the subject 




from some of the approximations that have usually been made. The 
latest and most distinguished writers upon this subject are General Char- 
bonnier and his pupil Sugot, both members of the celebrated Commission 
de Gavre (the latter its chief engineer) from which emanated all the knowl- 
edge of ballistics that was used during the war, not only for the French 
army but for ours. 
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The new methods here described consist, first, in the use of new experi- 
mental methods for obtaining the pressure in the gun, the position of the 
shot, and the muzzle velocity; second, in making the theory depend upon 
a new differential equation, and third, in the use of the graphic methods 
of calculation and integration that are common in engineering problems. 
The experimental methods hitherto have depended on the crusher gauge 
of Sir Andrew Noble and the Le Boulenge chronograph, _both almost 
unchanged since their invention fifty years ago, and both giving but a 
single value of the pressure or the velocity. To be sure the Sebert veloci- 
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meter gives a continuous record of the position of the shot as a function of 
the time. It will readily be seen how a continuous record, with an infinity 
of values of the pressure, as well as a more direct measure of the positions 
of the shot, will increase our knowledge of the whole process. 

We must first decide upon the equation of state that is to be used for 
the gases generated by the explosion. It has been shown in a previous 
paper in these Proceedings (Vol. 5, pp. 286-288, July 1919) what equations 
of state are compatible with the constancy of the difference of specific 
heats at constant pressure and constant volume, and we shall use the 
equation of Clausius, 
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\P + 



(v - a) = RT, 



(1) 



T(v + 0) 2 J 

in the simplified form, suitable for the high temperatures concerned (2000 °- 
3000° C) 

p(v - V )=RT (2) 

We shall use nearly the notation of Charbonnier and Sugot (Cours de 
Balistique, Imprimerie Nationale), but after having striven desperately 
with the multitude of engineers' units, pounds per square inch, kilograms 
per square centimeter or decimeter, inches or centimeters for bore and 
feet or meters for travel, and unknown units for velocity, acceleration and 
mass, we shall go back to first principle and express everything in C. G. S. 
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absolute units. It should be an axiom that no formula should contain 
any (arabic) numeral quantities except natural constants of the universe 
such as 7r, and that every formula should be correct in whatever system 
of units it is interpreted. It is hard to get engineers to agree to this. 
Consequently our unit of pressure shall be the dyne per square centimeter, 
with the convenience that a megadyne is between, and nearly equal to, 
the standard atmosphere and the kilogram per square centimeter at a 
place usually not specified, that v is the specific volume, in our units in 
cubic centimeters per gram or liters per kilogram, while T is the tempera- 
ture on the absolute scale, in our units centigrade degrees. If we con- 
sider a load of powder of mass w contained in a volume V, we define 
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the density of loading as A = w/V = 1/v, while if the fraction of 
the load z has been converted into gas we have v — V/coz. Accordingly 



RT 

v — V\ 



RT 



RTz 



V/z 



V 



zr, 



(») 



It is frequently assumed that during the burning the temperature of the 
gas remains constant, although this can hardly be the case during the 
expansion in the gun. Although the specific heats of the gases vary 
largely with the temperature, as shown by the classical experiments of 
Mallard and lye Chatelier (who found them to be linear functions of 
the temperature) we shall here consider them as constant, taking mean 
values. This hypothesis can be corrected later. 




If the powder is contained in a variable volume c and coz has been 
burned, while &>(1 — z) remains in solid form, of density 5, the gases will 
have the space V = { c — o>(l — z)/8 } to occupy, so that the specific volume 



will be v = 



- (1 — z) } /wz and we have 
8 



P 



RTuz 



(4) 



c — co(l — z)/8 — wzr\ 

as the relation between p, z, c, and T. Let 5 be the area of the bore, 
c' the volume of the chamber, and x the position or "travel" of the shot, 
so that 

c = c' + Sx. (5) 
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(6) 



In order to write the dynamical equation of motion of the shot we shall 
assume that the resistance of the rifling (sometimes known as the passive 
resistance) which is due to the friction of the shot against the lands, and 
would naturally contain a term proportional to the pressure, is a linear 
function of the pressure 

R = Ro + bp. (7) 




If then m is the mass of the shot, u its velocity, t the time, and if we take 
account of the inertia of the gas already formed, which is really uz, by 
means of a factor X, we have 

~ = u, (8) (m + Xco) ^ = Sip - p ) -Ro- bp, (9) 
at at 

where p Q is the atmospheric pressure, which is so small that it may gener- 
ally be neglected. 

By integration of (9) with respect to the time, 



(w + \co)u = (S - b)fpdt - iSp + R )t, 



(10) 
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(m + Xw)* = (S - b)fdtfpdt - (Sp + i? )* 2 /2. (11) 

o o 

If p is given in the form of a curve by the indicator, and the two integrals 
are evaluated for the whole time the shot is in the barrel, the determina- 
tion of the muzzle velocity and of the whole length of the barrel, every- 
thing else being known, will furnish two linear equations for the deter- 
mination of R and b. On the other hand if the times at which the shot 
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passes a number of points are known, or if x is experimentally given as a 
function of t either of the equations (10) or (11) gives a means of verifying 
the assumption (7) and as many points as may be for determining the 
constants by least squares. 

It will now be necessary to take account of the first law of thermo- 
dynamics. If U is the energy of the hot gas per unit of mass, dQ the 
amount of heat taken in, and Rdx the work done against friction, which 
is converted into heat, which stays in the gas, we have 

dQ + Rdx = d(azU) + pdV. (12) 
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If we now add (9) multiplied by (8) to (12) we have 

dQ = (m + \o))udu + d(cozU) + - pdz + Spzdx (13) 

and if there are H units of heat evolved for each unit of mass of powder 
burned, and dq lost to the gun by conduction, we obtain the equation of 
conservation of energy, 

Hwdz — dq — (m + \u>)udu + d(o>zlI) + - pdz -\- Sp dx (14) 

5 

Now if the specific heats are constant, we have 

U = CJ = -*I- = g£ ~ *> (15) 

K — 1 K — 1 

where k is the ratio of the specific heats. The quantity of heat Hdz — 
pdz/8 may be determined by calorimetric experiments in a bomb at con- 
stant volume, as is seen by putting dx and du equal to zero. Putting this 
equal to fdz/( k— 1) as is customary, neglecting q and p ot which can be 
done, and integrating (14) we have 

fwz , . M 2 wzRT , , 

- (m + Xw) + — (16) 



K ~ 1 ^ (C - 1 

which is Risal's equation, given in 1864. Combining this with (6) we 
eliminate the unknown temperature, which has never been experiment- 
ally determined, and obtain 

fd)Z — - (k — l)(m + Xco)m 2 . 
p = 2 ( 17 ) 

c' 4- S* — co(l — z)/5 — W772 

From this, if ^?, x, and m are experimentally determined, we may calculate 
z for each position of the shot. 

If we could neglect u 2 , or if it were proportional to z and if x, p, z, Were 
taken as coordinates, this would be the equation of an hyperbolic para- 
boloid, of which the sections z = const, are called isopyric lines by Lieut. 
Col. Hadcock in his paper on Internal Ballistics (Proceedings Royal Society, 
July 1, 1918 (479-509)). Col. Hadcock states that the expansion is 
neither adiabatic nor isothermal, but something between the two. This 
is true enough, as will be shown later, but it can hardly be justifiable to 
put, with any value of e, Col. Hadcock's equation (5) p(v — a) e /z e — K, 
for this makes the pressure proportional to z', instead of a linear fractional 
function of z, The assumption that u 2 /z is constant is not exactly true, 
but might be adopted without great error, as we shall see. 
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If we now put for brevity, 

-•- - - B, - = C, v - - = D, S - b = E, (18). 

co 5 <° 5 m + Xco 



Ro + 5^ 
S - b 



'•(^X^H 



where B and C are constants depending on the load, D a constant of 
the powder, E a constant of the gun and shot, and G a mixed mass and 
specific heat constant, or g = / — Gu 2 /z an approximate constant, we 
may write 

fz - Gx 2 



B + Cx - Dz 
B + Cx - Dz 



(19) 
(20) 



We have heretofore said nothing about the rate of burning of the powder. 
After the powder is all consumed, z = 1, and the expansion is adiabatic. 
We are interested in the period of combustion. It is now customary to 
assume, after Charbonnier, that the rapidity of burning of the powder 
is a function of the amount already burned, the function <p(z) being known 
as the form function of combustion, while a factor A is called the vivacity, 
being, for a given powder, inversely proportional to the linear dimensions 
of the grain of given form. The rate of burning is also proportional to 
some increasing function of the presssure, let us say P(p). The simplest 
functions, except linear functions, are some powers of the variables, not 
in this case integral powers. If the power is progressive we take <p(z) = 
3", if degressive (as in our powders), <p(z) = (1 — 2)". We will also put 
P = p a , where a is positive, as the burning is faster the greater the 
pressure. We accordingly put 

~ = A<p(z)P( P ) = A(l - zYp". (21) 

at 

We thus have the three simultaneous equations, 

dt= ^— = g? = du . (22) 

A<p(z)P(p) u E{p - F) 

It is obvious that the simplest choice for independent variable is z, accord- 
ingly we take 

dx = u (23 s du = E (p- F) , 

dz A<p(z)P(p)' dz A<p(z)P{p)" 

and eliminating u by differentiation, 
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dz 2 



1 du 
A<pP dz 



A, 



P<p' + V P' 



dx) 
E{ P - F) 



dx 

dz 



(11 H + A 

\P dz <p/' 



Now by logarithmic differentiation of (1.9) we obtain 

/ - 2Gu ----- C - - D 



1 dp _ 
p dz 
or by means of (23), (24) 



2Gu^ 
dz 



c dx 

dz 



(25) 



(26) 



fz — Gu 2 



B + Cx - Dz 



dz 



f - 2GE( P - F) 



dx 
dz 



C 



dx 
dz 



D 



fz - GA V-P 2 r 



/dxV 



B + Cx 



Dz 



(27) 



so that we have finally 

dH _ E(p - F) ■ _ dx 
dz 2 AVP 2 dz 






+ P 



P 



f - 2 GE(p - E) 



dx 
dz 



C 



dx 
dz 



D 



fz - GA 2 <p 2 P 2 \ 



(dxV 
\dz) 



B + Cx 



Dz 



(28) 



This is a sufficiently complicated equation of the second order and of 
the third degree in dx/dz, but it is exact, -and no simplifying assumptions 
have been made. If we know <p, P, p being obtained from equation 
(19), it gives x, the travel of the shot, in terms of z, the fraction of the 
powder burned. The equation (28) may be integrated graphically. When 
x is known as a function of z, p may be so determined, and then t and u 
from equations (22). Thus x, p, u, and z may all be found in terms of 
t, and the direct problem is solved. The equation (28) will be consider- 
ably simplified if we neglect G, which is equivalent to assuming the 
temperature constant, which is usually done. 

We shall however use the differential equation (28) or rather (25) to 
solve the inverse problem. If x, p, and u are experimentally given, from 
which z is calculated by (19), everything in (28) is known except P and <p. 
If we begin with an approximate value of <p, say (1 — z) & and an assumed 
value of p", equation (28) will be a linear equation for a and /?, from which 
by a few trials, the exact values may be obtained. 

We now come to the experimental portion of the work. The pressures 
in the rifle were observed by means of the gauge described in these Pro- 
cbSdings for July, 1919, the film being afterwards placed in the lantern 
and an enlarged tracing being made with a pencil on squared paper. From 
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this the calculations were made by my assistant, Dr. E- A. Harrington, 
who has carried out all the experimental work except the determination 
of the travel of the shot. In the figures (p. 649), curve 1 represents p,t for 
the shot used, while 2 is a curve of a reduced charge of 7/10 normal. The 
curves were then graphically integrated twice by counting the squares. 
Curves 3 and 4 show the first and second integrals J* pdt and fdtfpdt, 
respectively. The travel of the shot was determined by observing the 
time of contact of the bullet with the end of a coil of fine wire forced 
down the barrel to a measured point, an oscillograph of high frequency 
being used, and the time being observed on one of the revolving drums 
previously described. This was done about a year ago by my then assis- 
tant, Mr. H. C. Parker, and the result, shown in curve 5, is the average 
of a good many shots. The exact agreement of the two curves 4, 5 shows 
the propriety of the assumption about the resistance, and determines b 
as (S—b)/S = 0.894. In this case R = 0. We can also neglect p 
so that F = 0. 

The muzzle velocity and the time of the shot leaving the barrel were 
determined by the method described in these Proceedings, April, 1920, 
developed in this laboratory by Thompson, Hickman, and Riffolt. The 
constants of the gun were as follows : 

Volume of chamber, c = 4.18 cm. 
Area of bore, S = . 456 cm. 

As for the constants of the powder, the density was determined by Dr. 
Harrington in the pyknometer as 1.58 gm./cm. No attempt was made 
to determine -q experimentally, that being best done in the bomb, but the 
value jj = 0.95 was assumed. The load was co = 3.100 gm. To deter- 
mine/, the "force" of the powder (which has the dimensions of energy/ 
mass) two methods were used, which agreed very well. If it is assumed 
that the point of inflexion on the observed pressure curve marks the end 
of combustion, z = 1, / is determined from equation (19). If, on the other 
hand, equation (19) is used with observed values of p , x, and u, it will 
give values of z increasing to a maximum, and then diminishing. This 
maximum should be equal to unity, and after a few trials a value of / 
can be found that makes it so. As stated, this method agrees with the 
other. We obtain / = 1.087 X 10 10 cm. 2 /sec. 2 , corresponding very well 
with the values for the French powders and cordite. (The values given 
by Charbonnier and Sugot involve gravitation units.) We thus obtain 

B = 0.716 cm. 3 /gm., C = 0.147 cm. 2 /gm., 

D = 0.32 cm.Vgm., E = 0.394 cm. 2 /gm. 

If we take k = 1.2, and X = 0.2, we have G = 0.333. Also we find 
that g (eq. 20) varies by about 16%, nearly as a linear function 
of x or z. 
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In order to show the functional relation (19) between the five quantities 
p, x, u, z, G, an "abaque a glissement" was constructed, so that by placing 
a transparent index on four of the variables the fifth could be found. 
This will be described elsewhere, but inasmuch as it was found to take 
longer to set the index than to set up the numbers on the arithmometer, 
of which we have several models, it was not used. After z was calculated, 
curves were drawn : 6(p,x), the indicator diagram (p. 652) for the gun; 7 (2, i), 
8(dz/dt,t), 9(x,z), 10(dx/ds,z), 11(1— 2)" for various values of n; 12(p,z), 
lZ{dp/dz,z) (with a different zero of ordinates); 14(#,z); 15(du/dz,z); 
16(dz/dt,z). 

For graphical differentiation, a simple contrivance was used consisting 
of two metal rulers jointed together at one end. Near the other end 
of one was placed a plane mirror normal to the edge of the rule. By 
looking at the curve in this, and making its reflection join it without a 
corner, the direction of tangency could be found more accurately than 
by making the ruler tangent. The other ruler is then set parallel to the 
horizontal lines of the paper, and for a horizontal distance of ten or twenty 
centimeters the vertical distance between the rules read off, giving the 
trigonometrical tangent. Thus in a few minutes we may lay down the 
derivative curve with nearly three figure accuracy. 

I have the following criticisms to make on the curves. After they had 
been given to the photographer it was discovered that the left-hand 
portions of (9) and (10) are incorrect, since dx/dz must vanish when x 
and u are zero. This is probably due to an error in choosing this point, 
which should be at the inflexion of the p curve, which should have been 
determined with the mirror directly on the film, instead of on the enlarge- 
ment. When this is done, the pressure to start the shot is found to be the 
same as on curve 2, which (like 1) is so good that films of two different shots 
may be laid one over the other without being distinguishable. Making 
this correction will slightly displace curves (3) and (4). I realize that I 
have given hardly enough here to verify the theory, but as this paper is 
too long already, and as I am very shorthanded for assistance, I consider 
it best to publish it for the method, and to leave the verification, which 
should be by the construction of the observed curves from calculation, 
for another paper. By the courtesy of Rear-Admiral Ralph Earle, we 
have a Naval one-pounder gun, to which we hope to apply the methods, 
on which we hope to improve by experience. For a larger gun, where 
the wires might be blown out of the gun by the blast, we have another 
method of finding the velocity in the bore. 

It should have been stated in our former paper, and is now stated 
here, that these experiments have been made possible by means of a 
generous grant from the Rumford Fund of the American Academy of 
Arts and Sciences, and to my colleagues on the Rumford Committee 
I hereby express my hearty thanks. 



